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ABSTRACT 

Impulsive nanoflares are expected to transiently heat the plasma confined in 
coronal loops to temperatures of the order of 10 MK. Such hot plasma is hardly 
detected in quiet and active regions, outside flares. During rapid and short heat 
pulses in rarified loops the plasma can be highly out of equilibrium of ionization. 
Here we investigate the effects of the non-equilibrium of ionization (NEI) on the 
detection of hot plasma in coronal loops. Time- dependent loop hydrodynamic 
simulations are specifically devoted to this task, including saturated thermal 
conduction, and coupled to the detailed solution of the equations of ionization 
rate for several abundant elements. In our simulations, initially cool and rarified 
magnetic flux tubes are heated to 10 MK by nanoflares deposited either at the 
footpoints or at the loop apex. We test for different pulse durations, and find 
that, due to NEI effects, the loop plasma may never be detected at temperatures 
above ~ 5 MK for heat pulses shorter than about 1 min. We discuss some 
implications in the framework of multi-stranded nanoflare-heated coronal loops. 

Subject headings: Sun: corona - Sun: X-rays 



1. Introduction 



Nanoflares - small scale highly transient h eating episod e s - are among the main can 



didates as source of coronal loop heating (e.g., iParker 



19881 : ICargilllll994l : iKlimchukl 12006 
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and references therein). The conjecture is still under debate because nanoflares have been 
hardly detected so far. There are many possible reasons for this difficult detection. For 
instance, a very frequent occurrence may inhibit the resolution of the single event. Also 
the efficient thermal conduction and low emission measure during the pulses can reduce and 
delay the signatures of the heating (e.g., iPeres et al.lll987l ; iReale fc Pereslll995l ). Also very 
small pulses distributed in a very finely structured loop may be difficult to detect. One of the 
main arguments invoked as crucial evidence of nanoflare heating is the d etection of high tem- 
perature (10 MK) plas ma components in observations of coronal loops (jCargill fc Klimchuk 
19971 ; iKlimchukl 120061 ) . outside of proper flares. The presence of such hot component is 
often predicted by hydrodynamic modeling of coronal loops heated by transient pulse s 
dPatsourakos & Klimchukl I2OO5I : IPatsourakos & Klimchukl bood : ICargill & Klimchukl I2004T ) . 
In particular any heat spike able to bring the loop to the observed brightness should be so 
intense as to heat the plasma to temperatures of the order of 10 MK at least for a transient 
time interval. In the hypothesis of a finely structured loop where a whole distribution of heat 
pulses occur continuously, such plasma might be detectable, at least as a hot tail in the emis- 
sion measure distribution. The evidence of hot plasma has been difficult so far. For instance, 
a dominant plasma component at about 3 MK is shown by recent thermal maps of active 
regions obtained from wide- band multi-filter i maging observations (jReale et al.l 120071) with 



the X -Ray Telescope (XRT, iGolub et al.l 120071 ) on board the Hinode mission ( IKosugi et al. 
2007h . 



There are possible explanations also for the difficult detection of hot plasma. The most 
immediate one is linked to the inertia of the plasma dynamics. A heat pulse deposited in 
a coronal loop drives evaporation of chromospheric plasma. The loop is filled with hot and 
dense plasma which makes it bright in the X-rays. For short heat pulses, the plasma may 
evaporate from the chromosphere on time scales longer than the heat duration. Therefore 
the heated strand may only later acquire enough emission measure to become visible, while 
it is cooling. Ho wever, in sho rt loops (~ 10 9 cm) the sound crossing time is of the order of 



one minute (e.g., iReald 120071 ) and we may expect to detect hot spots within less than half a 
minute. 

There is another less obvious effect which may make the detection of hot plasma harder: 
the time lag of the plasma to change its ionization from a cool to a hot state. An impulsive 
energy input drives plasma thermal and dynamic changes on relatively short timescales. 
Electron excitation, de-excitation, ionization and recombination processes of the ion species 
have other timescales. If the timescale, for instance, of the temperature evolution is much 
shorter than the ionization and recombination timescales, the degree of ionization can be 
yery different from the equilibrium conditions corresponding t o the local electron temperature 
(jShapiro &: Mooreill977l ; lOrlando et al.lll999l ; IKlimchukl 120061 ) . So, during a fast temperature 
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increase, the plasma ions can be at a lower ionization state than the equilibrium state 
corresponding to the instantaneous temperature. Such non-equilibrium of ionization (NEI) 
effects may become important in the interpretation of what we observe. If the heat pulse 
is sufficiently short, the ions and their emitted radiation may not even have enough time to 
"sense" the hot temperature status, they would adjust to the temperature variations deep 
in the cooling phase, and we would detect radiation from cooler ion conditions at any time. 



NEI effects in nanoflare-heated loops have already been investigated in the past. iMariska et al 



(119821 ) found that NEI can significantly alter the relative ionic abundances in the quiet 



Golub et al . (1989) e xami ne the effect of NEI on the observability of coronal vari- 



Bradshaw fc Mason! (20031) found that NEI can modify considerably the radiative 



corona, 
at ions. 

losses function. iMuller et al.l (120031 ) inve stigated the role of NEI in the transition region 
brightenings driven by loop condensations. iBradshaw &: Cargilll (120061 ) model nanoflare heat- 
ing in coronal loops and remark the importance of NEI effects, the presence of hot plasma 
with low emission measure, and Doppler-shifts as possible diagnostics of the heating. 

Here we investigate the effect of NEI during nanoflaring activity on the detectability 
of hot plasma in coronal loops. We will model a coronal loop strand heated by nanoflares 
of different durations, compute the corresponding evolution of the relevant ion species, and 
compare it with the evolution in full ionization equilibrium. We will evaluate the effects 
on the expected temperature distribution of the loop emission measure, which will tell us 
about the existence of significantly emitting hot plasma components. This will allow us to 
put constraints on the nanoflare characteristics which lead to the observability of hot plasma 
and indirectly also on the fine loop structuring. In Secj2]the modeling is described, in Secj3] 
the results are illustrated and discussed in SecJH 



2. Modeling 

Our model is set up to explore the conditions for plasma confined in a coronal magnetic 
flux tube and heated to 10 MK by nanoflares to be detectable as hot plasma considering 
the effects of NEI. The coronal loop contains low f3 plasma, and, as customary for standard 
loop models, we assume that plasma moves and transports energy only along the ma gnetic 



field lines, so that a one-dimensional description is adequate (e.g., iPeres et al.lll982l ). The 
model takes into account the gravity stratification, the thermal conduction (including the 
effects of heat flux saturation), the radiative losses, an external heating input, and the NEI 
effects. The following fluid equations of mass, momentum, and energy conservation are 
solved, considering only the relevant components along the loop magnetic field lines: 
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dp 
dt 



+ V • pv = 



<9pv 
~dt 



V ■ pvv + VP = pg 



(2) 



<9p£ 
pv 



+ V • (pE + P)v = 

■g-V-g + Q(s,t)-n e n H A(T) 



(3) 



where E 



is the total gas energy (internal energy, e, and kinetic energy), t is the time, s is the coordinate 
along the loop, p = pm^nn is the mass density, p = 1.26 is the mean atomic mass (assuming 
solar abundances), m# is the mass of the hydrogen atom, «h is the hydrogen number density, 
n e is the electron number density, v the plasma flow speed, P the pressure, g the gravity, 
T the temperature, q the conductive flux, Q(s,t) a function descr ibing the transient input 
heating, A(T) is the radiative losses per un it emission measure (e.g. iRaymond fc Smithllf977l : 



Mewe et al. 1985; Kaastra Mewe 



20001 ). Here we consider the radiative losses function in 



equilibrium of ionization. Although this assumption makes our description not entirely self- 
consistent, it does not affect significantly our results, because in the critical evolution phases 
the energy losses are dominated by transport by thermal conduction (see Sec J3.ll) . We use 
the ideal gas law, P = (7 — l)pe, where 7 = 5/3 is the ratio of specific heats. 

The set of the continuity equations for each ion species is expressed as: 



-^ + V-nfv = Rf 



Nelem) 



1 N Z ) 



(4) 



where Rf 



n e [n i+1 a i+1 



+ nf_ 1 Sf_ 1 -nf(af + Sf)\ 



nf is the number density of the i-th ion of the element Z, N e i em is the number of elements, 
Nf on the number of ionization states of element Z, af are the collisio nal and dielectr onic 



recombination coefficients, and Sf the collisional ionization coefficients (I Summers! 1 1 9 741 ) 
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Given the importance of rapid transients, fast dynamics and steep thermal gradients 
in this work, we consider both the classical and saturated conduct i on reg imes. To allow 
for a smooth transition between them, we follow iDalton &: Balbud (119931 ) and define the 
conductive flux as 



1 1 

Qspi Qsat 



(5) 



Here q sp i represents the classical conductive flux (jSpitzerlll962l ) 



<?s P i = -k(T)VT 



(6) 



where the t hermal conductivity is k(T) = 9.2 x 10 7 T 5//2 erg s 1 K 1 cm 1 . The saturated 
flux, g sat , is (jCowie fc McKedll977l ) 



g sa t = -sign (VT) 50pc; 



(7) 



where c s is the isothermal sound speed, and is a correction factor of the order of unity . 
We set (ft = 1 according to the values su ggested for the coronal plasma (IGiulianil Il984l ; 



Borkowski et al.l Il989l . iFadeyev et al.l |2002| . and references therein) . 



The transient input heating is described empirically as a separate function of space and 



time: 



Q{s,t) = H x g(s) x f(t) 



where g(s) is a Gaussian function: 



g(s) = exp[-(s-soY/2a 2 ] 



(9) 



and f(t) is a pulse function: 



0, t < 



f(t) = { 1, < t < t H 



(10) 



0, t>t H 
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The calc ulations described in this paper were performed using the 1-D version of the 
FLASH code (IFryxell et al.ll2000l ). an adaptive mesh refinement multiphysics code. For the 
present application, the code has been ext ended by additiona l computational modules to 
handle the plasma thermal conduction (see lOrlando et al.ll2005l for details of the implemen- 
tation), the NEI effects, the radiative losses, and the heating function. The implementation 
of the description of the ionization balance in the flash code is described in the Appendix. 

We consider an initially cool and rarified semicircular magnetic flux tube (loop or loop 
strand) with half-length L = 3x 10 9 cm, and uniform cross-section area. The area appears 
only as a multiplicative factor in this description and the size is typi cal of active region 
loops. The loop is initially at e quilibrium accord ing to loop scaling laws fjRosner et al.lll978l ) 
and to hydrostatic conditions flSerio et al.lll98ll ) assuming it lies on a plane vertical to the 
solar surface. The initial base pressure is p$ = 0.055 dyn/cm 2 , corresponding to a maximum 
temperature To = 0.79 MK, at the loop apex. The corona is linked with a steep transition 
region to an isothermal chromosphere uniformly at T c = 20000 K and 0.5 x 10 9 cm thick. 
We assume that the loop is symmetric with respect to the vertical axis across the apex, and 
therefore we simulate only half of the loop, with a total extension of 3.5 x 10 9 cm. All the ion 
species are assumed initially in equilibrium of ionization in the whole computational domain. 

At the coarsest re solution, the adaptive mesh algorithm used in the flash code (paramesh; 
MacNeice et al.ll2000l ) uniformly covers the computational domain with a mesh of 16 blocks, 
each with 8 cells. We allow for 3 levels of refinement, with resoluti on increasing twice at 
each refinement level. The refinement criterion adopted (jLohnerlll987l ) follows the changes in 
density and temperature. This grid configuration yields an effective resolution of ~ 3.4 x 10 6 
cm at the finest level, corresponding to an equivalent uniform mesh of 1024 grid points. 
We use fixed boundary conditions at s = and reflecting boundary conditions at s = s max 
(consistent with the adopted symmetry). 



3. Results 

3.1. The Simulations 

In the loop outlined above we inject one nanoflare as intense as to heat the plasma 
and keep it at ~ 10 MK, if the heating were steady. The choice of the parameters is 
dictated by our scope of exploring the influence of NEI on the detectable thermal conditions. 
In this perspective, the pulse duration becomes the critical parameter: if the heat pulse 
is short enough, the plasma may not have enough time to adjust to ionization conditions 
appropriate of 10 MK, before the heating is off. Therefore, we consider longer and longer 
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heat pulse durations with a logarithmic sampling, i.e. tu — 5 s, 30 s, 180 s. For each of 
these durations we consider two possible locations of the pulse depositions: at the footpoints, 
namely Sq = 8 x 10 8 cm from the base of the chromosphere (i.e. 3 x 10 8 cm from the base 
of the corona) and at the apex, namely So = 3.5 x 10 9 cm. The spatial width of the pulse 
is smaller at the footpoints (er = 10 8 cm) and larger at the apex (er = 5 x 10 s cm). To have 
the same total energy deposition rate, the maximum volume deposition rates are Hq = 1.5 
erg cm -3 s _1 and H = 0.2 erg cm -3 s _1 , respectively. 

We come out with a total of six simulations. They are carried out for a time which 
covers the pulse duration and a couple of loop cooling times: 



Tde 



120- 



360 s 



where L g = 3 is L in units o f 10 9 cm and T7 m 1 is the loop maximum temperature in units 
of 10 MK (ISerio et al.lll99ll ). All simulations then span a time interval of about 800 s. 



The e yolution of nanoflaring plasma confined in coronal loops is well-known from previ- 



ous work ( 


Peres et al. 


1993 


Warren et al. 


2002' 


Warren et al. 2003:Patsourakos & Klimchuk 


2005 Testa et al. 2005) 


and is on a smaller scale similar to that of properly flaring loops (e.g. 


Naga|l980 


Peres et al.lll982 


Doschek et al.lll983 


Naeai & Emslie 1984: Fisher et al. 1985: Peres et al. 


1987l]Reale & Peres! 1995h . As additional feature, our simulations include the effect of satu- 



rated thermal conduction, which might be i mpor tant here since we study the evolution on 



small time scales (see also iKlimchuk et al.l 120031 ). In the following, we just draw a basic 



outline of the results adapted to our simulations. 

Let's consider the simulation with tn — 5 s and heat pulses deposited at the footpoints 
(Fig. [1]). The sudden heat deposition determines a local increase of temperature (to ~ 10 
MK) and pressure (to ~ 1 dyn/cm 2 ). A fast thermal front propagates upwards along the 
loop. The cool chromosphere is heated and expands upwards with a strong evaporation 
front. By the time the thermal front has reached s ~ 1.5 x 10 9 cm and the evaporation 
front s ~ 6 x 10 8 cm, the heat pulse is already over. Due to efficient thermal conduction, 
the plasma then immediately begins to cool, already during the propagation of the thermal 
front, and the maximum temperature rapidly decreases to ~ 2.5 MK in less than a minute. 
The thermal front reaches the apex, and therefore the loop thermalizes, in ~ 30 s. The 
impulsive evaporation front moves at velocity of about 300 km/s (FigJHO) and fills the loop 
in ~ 75 s. The density keeps on increasing throughout the loop with more gentle fronts (the 
peak velocity rapidly decreases to less than 200 km/s) for further tens of seconds. Meanwhile 
the plasma accumulates at the apex reaching a density of about 10 9 cm -3 . The compression 
heats the plasma again above 3 MK at the apex. Then the plasma begins to drain and the 
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Fig. 1. — Evolution of the plasma temperature (A), density (B), pressure (C), and velocity 
(D) distributions along the loop from the chromosphere to the coronal loop apex, sampled 
every 5 s from to 100 s (thinner and thinner lines), for heat pulse location at the loop base 
and heat pulse durations tjj = 5 s. 



density to decrea se (not shown), following the radiation cooling time (jCargill fc Klimchuk 



2004 : iReald 120071 ) 



The overall evolution does not change much in the other simulations. Longer heat 
pulses drive more and longer plasma evaporation. Pulses deposited at the apex produce 
downward thermal fronts but as soon as the thermal front hits the chromosphere, it drives 
an evaporation front similar to that driven by the pulses at the footpoints. 







10 20 
Ionization State 



30 



Fig. 2. — Ionization/recombination timescale (in log scale) for various ion species (dots) 
vs the ionization state for the labeled elements, computed for a temperature 10 MK and a 
density 10 8 cm -3 . The dashed horizontal line marks 100 s. 



3.2. Non-equilibrium of Ionization 



The importance of highly transient processes such as NEI is basically dictated by their 
timescales related to the timescales of the dynamics and heating/cooling driven by the 
nanoflares. Fig. [2] shows the combined ionization/recombination timescale for the i-th ion 
species of the element Z, derived from Eq.(j4]) as 



T i ~ { Z i oZ\ ' ("^) 

n e {af + Sf) 



for twelve important elements, computed for a temperature 10 MK and a den s ity 10 8 cm 
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Analogous timescales have been provided, in tabular form, by iGolub et al.l (119891 ) in the 
study of observable variability of spectral lines in soft X-ray and XUV regions of the solar 
corona. 

Clearly, r-f increases with the ion species and becomes larger than 100 s for the highest 
ion species of C, N, O, Ne, Mg, Si, S, Ar and Ca. More species are involved in high 
ionization times for Fe and Ni. The modeling shows that the plasma immediately cools 
down considerably as soon as the heat pulse stops. As a consequence, from Fig. [2] we expect 
that the high ion species will have no time to adjust to high temperature status if the heat 
pulse lasts significantly less than 100 s. 



The ionization and recombination timescale is dominated by either ionization or re- 
combination rate, depending on the sign of the temperature jump and on the ionization 
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state of the element considered: in general, is dominated by the ionization rate for 
the lowest ion species and by the recombination rate for the highest ion species. The ion- 
ization state for which ionization and recombination rates are comparable depends on the 
temperature: the higher the temperature, the higher is this ionization state. In the case 
cons idered here, is dominated by the ionization rate for most of the ion populations (see 



also iPatsourakos fe Klimchukl 120061 ) and by recombination rate for the top ionization states 
(e.g. Si XV, S XVI, Ca 20, Fe XXVI). As shown later, we found the effective ionization tem- 
perature generally lower than the electron temperature, even when the plasma is cooling; 
as a consequence, in the case of nanoflares discussed here, is generally dominated by the 
ionization rate for most of the ion populations. 

As an example, Fig. [3] shows the distribution of population fractions of Fe along the 
loop derived assuming equilibrium ionization (upper panels) and considering the deviations 
from equilibrium ionization (lower panels) for the simulation with heat pulses deposited at 
the footpoints lasting tg = 5 s (see also Fig. [TJ. During the early phase of the evolution, 
the deviations from equilibrium of ionization are very large at the loop footpoints due to 
the sudden local increase of temperature to ~ 10 MK (see left panels of Fig. [3]). Whereas 
significant populations are expected above Fe XX low in the loop (below 1.5 x 10 9 cm, reddish 
lines in the upper left panel), no such highly ionized Fe species appears when NEI is taken into 
account. As expected, the high ion species are not able to adjust to the plasma temperature 
before the heat pulse is over. Then the plasma cools down reducing the deviations from 
equilibrium ionization. The deviations are still significant later at t = 30 s, when the loop 
has thermalized (see Fig. [TJ and the maximum temperature is ~ 3 MK. At this time, no 
highly ionized Fe species is present both considering and not considering NEI, because of 
the effective cooling. 

This time lag of the evolution of the ion species has important implications on the tem- 
perature diagnostics of the plasma from its emission. To address this issue, from the modeled 
population fractions, we derive the temperature which best matches the actual ionization 
state of the ion species for each of the elements in our simulations: operatively this is found 
as the one of the equilibrium temperatures with the most similar three most populated ion- 
ized states. It is worth noting that at least two population fractions are necessary to find a 
unique value of temperature and that this "NEI" temperature should not be intended as an 
exact temperature but as the temperature which best describes the ionization state. This is 
done at each time and position along the loop. Fig. 2] shows the evolution (on a logarithmic 
time scale) of the plasma maximum temperature and of the maximum NEI temperature for 
five representative elements (C, O, Mg, S, Fe). In the same figure the evolution of the coronal 
emission measure at the maximum plasma temperature is also shown for comparison. The 
ion species adjust to a hotter status very gradually, on a time scale of about 100 s (Fig. [2]). 
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1.000 F 




Fig. 3. — Distributions of Fe population fractions (color-coded on the right axis of the upper 
right panel) along the loop from the chromosphere to the coronal loop apex, assuming equi- 
librium ionization (upper panels) or considering the deviations from equilibrium ionization 
(lower panels), at the labeled times (5 s and 30 s), for the simulation with heat pulses located 
at the loop base and heat pulse duration tu = 5 s. 

Since this time is slightly longer than the time taken by the emission measure to increase 
significantly (as shown by the bottom panels), we expect observable effects in the X-ray 
band. For i# = 5 s, despite the maximum electron temperature is above 10 MK, all the ion 
species are never represented by a temperature larger than about 3 MK all over the event, 
for any pulse location. For tu = 30 s, Fe (and Ni) reaches a temperature about 7 — 8 MK 
for pulses located at the footpoints and about 6 MK for pulse deposited at the apex. This 
occurs for a very limited time range (between 100 and 200 s). The other elements mostly 
stay at temperatures well below 5 MK. The results change quite substantially for tn = 180 
s. Although with delay (> 100 s), here many elements have enough time to adjust to ioniza- 
tion stages typical of higher temperatures, between 7 MK and 10 MK, for both heat pulse 
locations. The high ionization state is maintained for longer time intervals (more than 200 
s). 
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The effect of long-lasting NEI is illustrated synthetically in Figs. [5H3 which shows the 
total loop emission measure distribution with temperature, EM(T), averaged over different 
time intervals. The figure compares EM(T) obtained either taking or not taking the devi- 
ations from equilibrium ionization into account; the former is obtained using the average 
"NEI" temperature among those of all the elements considered, the latter using directly the 
electron temperature. 

Averaging over the heat pulse duration only, we can see that EM(T) with and without 
NEI are invariably very different. The shorter the pulse duration, the more the hot compo- 
nents are missing from the NEI EM(T) distributions. For tn = 5 the hottest components 
including NEI are only at about 1 MK, which increases to 3 MK for tu = 30 s. The dif- 
ference is less considerable for tn = 180 s, but even in this case the hot components are 
hardly above 5 MK, against 10 MK without NEI. The differences are unsurprisingly reduced 
when we average over longer times, which should approach a situation closer to the realistic 
observations. Over a few durations of the heat pulses the differences are still significant for 
the shortest heat pulses: for tn = 5 s, with NEI we do not find components above 2 MK, 
whereas without NEI there are components almost to 10 MK. For tu = 30 s the two EM(T) 
distributions are more similar, but with NEI we do not find components hotter than 6 MK. 
For the longest pulse duration, the EM(T) distributions almost coincide. 



4. Discussion and conclusions 



The target of this work is to explore the effect of non-equilibrium of ionization (NEI) 
in nanoflare-heated loops. The import ance of NEI e f fects on coronal observations had been 
already pointed out and evaluated by iGolub et al.l (119891 ) . That study was devoted more 
specifically to the diagnostics of rapid variability. Here, we focus on the effects on the 
detection of hot plasma in nano-flaring loops. To this purpose we set up hydrodynamic 
simulations of plasma confined in an active region loop heated to ~ 10 MK by more or less 
long nanoflares. As new achievements specifically set up for t his study, we have inc luded 
the effect of the saturated thermal conduction in the modeling (IKlimchuk et al.ll2003l ). and 
computed the evolution of the ion population fractions of several important elements driven 
by a short heat pulse. Each ion species takes a characteristic time to reach equilibrium 
conditions for a certain plasma temperature and density. If the plasma conditions vary 
on very small time scales, the ions are unable to adjust to the new and rapidly changing 
conditions and the emitted spectrum during the plasma variations can become very different 
from the one expected in equilibrium conditions. If the variation is a sudden heating followed 
by a sudden cooling, the emission may never or only partially adapt to the transient hot 
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conditions. To evaluate the importance of this effect and its dependence on the heat pulse 
parameters, we have computed the actual temperature of the plasma as it would appear in 
the radiation spectrum according to the effective ionization state. 

We find that the effects are significant if the heat pulses last less than a minute or so. 
For durations of a few seconds the spectra will never show temperatures higher than 2-3 
MK even though the electron temperature overcomes 10 MK. This effect is little dependent 
on the location of the heat pulse, except for details. These results can be considered valid 
in wide generality in spite of some limitations of our approach. Our assumption of pulses 
leading to 10 MK is, of cou rse, arbitrary, although this ha s been taken as possible typical 



temperature in other works (IPatsourakos fc Klimchukll2006l ). Conclusions are certainly valid 



for less intense pulses, which heat less the plasma. We do not expect dramatic differences 
also for more intense pulses (e.g. 20 MK), because the thermal conduction becomes even 
more efficient and the plasma initial cooling (after the end of the pulse) will be faster. Also 
the details of the shape of the heat pulse should not have much influence, provided that 
the overall duration is the same. In general, of course, more gentle heating should make 
plasma reach equilibrium conditions more easily, while the opposite occurs for more spiky 
pulses. More important is instead our assumption of a single pulse inside a given thread. 
With it, we are excluding a repetition of the pulse inside the same thread, or, at least, a 
very long repetition time, so long as to have a rarified loop (~ 10 8 cm -3 ) again before the 
new ignition. More frequent pulses would be released in a denser plasma which is faster to 
adjust to ionization equilibrium, and in which, therefore, we expect less NEI effects. The 
ionization times in fact scale inversely with plasma density. Heat pulses released in a ~ 10 9 
cm -3 dense plasma, i.e. about ten times denser than our initially rarified flux tube, would 
excite highly ionized element species about ten times more rapidly, with typical times of 
about ten seconds (Fig. [2]). 

In the range 1 — 10 MK the plasma emits radiation mostly in the X-ray band and the 
spectra are dominated by the emission lines produced by the reco mbination and deexcitation 



of highly ionized atoms (e.g., iTucker fc Korenlll97ll . lKatolll976l ). For instance, hundreds of 



lines of at least ten Fe ion species populate the spectra of plasma in that range of tem- 
perature. For this reason, the ionization status of the emitting ions becomes extremely 
important in the thermal "appearance" of the plasma, which we detect through its radiation 
with remote sensing, while the continuum electron bremsstrahlung emission - which adjusts 
immediately to thermal and dynamic plasma variations - becomes more important at much 
higher temperatures. Our results therefore apply both to the analysis of single lines through 
high resolution spectroscopy, and to the diagnostics from wideband multi-filter instruments 
in the X-ray and EUV bands, e.g. Hinode/XRT or the Atmospheric Imaging Assembly 
(AIA) on board the forthcoming Solar Dynamic Observatory. 
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According to these considerations, a hard detection of hot plasma may point to a scenario 
of coronal loops heated by single short and relatively intense nanoflares. Very short durations 
and very long repetition times naturally imply that the pulses should be deposited always 
in different strands and that, therefore, the strands involved should be a very large number, 
i.e. the loop must have a very fine transversal structure. A scenario of loops heated by 
short nanoflares has some important implications. Hot plasma may be difficult to detect not 
because of a limitation of the telescopes, detectors and filters, but due to an intrinsic property 
of the plasma emission. The critical parameter for detectability of hot plasma is the duration 
of the heat pulse and ~ 1 minute is the critical duration value, a useful indication for detailed 
models of nanoflaring mechanisms. The fine temporal structure would also automatically 
imply a very fine spatial structure, difficult to resolve. 

This picture is fully compatible with the current scenario of nanoflaring multi-stranded 
coron al loops. For instance, it naturally involves t he hot-underden se / cool-overdense loop 
cycle ( 1 Warren et al.l l2002|jCargill fc Klimchukl l2004j . iKlimchukl 120061 ) . Less trivial is the ef- 

fect on the diagnostics of the detailed thermal structure along and across the loops (e.g., 

Cargill fc Klimchukll20o4 kschwanden fc Nightingale!l2005l . Illeale fc Ciaravella!l2~0~0~6l . ISchmelz et~al 
20071 ). which we defer to later work. 



In conclusion, we have investigated in detail the effect of the non-equilibrium of ioniza- 
tion driven by nanoflares in finely structured coronal loops on the detection of hot plasma 
in the loops, and provided constraints on the heat pulse duration which may make or make 
not the detection difficult. Further more detailed information and diagnostics are expected 
from high resolution spectroscopy such as that from EUV Imaging Spectrometer on board 
the Hinode mission or from the Extreme Ultraviolet Variablity Experiment (EVE) on board 
the forthcoming Solar Dynamic Observatory. 
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A. Non-equilibrium ionization in the FLASH code 



The non-equilibrium ionization is added to the FLASH code using a method of time 
splitting between the hydrodynamic and the NEI numerical modules. A fractional step 
method is required to integrate the equations and in particular to decouple the NEI solver 
from the hydrodynamic solver. For each timestep, the homogeneous hydrodynamic transport 
equations given by Eqs. [THU are solved using the FLASH hydrodynamic solver with R = 
0. After each transport step, the stiff system of ordinary differential equations for the NEI 
problem with the form: 



dnf 
~dt 



N, 



spec) 



(Al) 



are integrated. This step incorporates the reactive source terms. Within each grid cell, the 
above equations can be solved separately with a standard ODE method. Since this system 
is stiff, it is solve d using t h e Bad er-Deuflhard time integration solveiQ with the MA28 sparse 
matrix package. iTimmed (119991 ) has shown that these two algorithms together provide the 
best balance of accuracy and overall efficiency. Note that the source term in the NEI module 
is adequate to solve the problem for optically thin plasma in the coronal approximation; just 
collisional ionization, auto-ionization, radiative recombination, and dielectronic recombina- 
tion are considered. 
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Fig. 4. — Evolution of the maximum electron temperature (solid lines) and of the maximum 
"NEI" temperature (see text for definition) for the labelled elements (dashed and dotted 
lines). The figure shows the results for heat pulse location at the loop base (left panels) and 
at the apex (right panels) and for heat pulse durations tn — 5 s, 30 s, and 180 s. Bottom 
panels: evolution of the coronal emission measure at the maximum plasma temperature. 
The time is in log scale to zoom up the fast initial evolution. 
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Fig. 5. — Distributions of emission measure vs temperature averaged over different time 
intervals obtained either taking (thick lines) or not taking (shaded areas) the deviations 
from equilibrium ionization into account (see text) for heat pulse location at the loop base 
(left) or at the apex (right) and heat pulse duration tu = 5 s. 



-21 - 




10 6 10 7 10 8 10 6 10 7 10' 

Temperature [K] Temperature [K] 



Fig. 6. — Same as Fig. [5l for heat pulse duration t H = 30 s. 
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Fig. 7. — Same as Fig. [5j for heat pulse duration t H = 180 s. 



